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SUMMARY 


A  method  for  producing  normal  random  variables  in  terms 
of  uniform  random  variables  U^,  If  Y  =  +  U^, 

then  choosing  one  of  the  four  random  variables  2Y  -  3,  (4T  -  6)/3 
(Y  -  7)/2  or  (Y  +  4)/2  in  the  proportions  •8635,  #11506, 

•00372  and  #00372  will  produce  the  reruired  normal  variate  98#6 
per  cent  of  the  time#  The  other  1#4  per  cent  is  devoted  to  the 
teiil  or  a  rejection  technique  in  order  that  the  composite  be 
exact#  The  method  leads  to  ver7  fast  computer  programs  v^ich 
are  easy  to  code  and  occupy  little  space  in  the  computer. 


THE  METHOD 


The  following  method  for  generating  normal  random  variables 
In  a  computer  Is  similar  to  that  reported  In  reference  [2J.  It 
is  slightly  faster  and  simpler*  Programs  based  on  this  method  take 
about  half  again  as  much  time  as  the  very  fastest  methods,  [l]* 

The  fastest  methods,  however,  require  a  few  thousand  storage 
locations  and  are  relatively  complicated*  The  basic  normaJ  random 
variable  generator  for  a  largo  computer  probably  should  use  the 
fastest  available  method,  but  in  situations  where  ease  of  pro¬ 
gramming  or  storage  requirements  are  more  important,  the  following 
method  is  particularly  well  sill  ted* 

We  want  to  generate  a  normal  random  variable  X,  density 
—  5  —  5x  ^ 

(2n)  *■'^6  *  ,  by  using  independent  random  variables  U^,  U^,**., 

uniformly  distributed  over  the  interval  (0,l)*  We  assume  our  com¬ 
puter  has  a  subroutine  for  providing  a  satisfactory  sequence  of 
U's*  The  idea  is  to  use  the  sum  of  3  tiniform  variables, 
y  =  +  U^,  in  one  of  the  four  linear  forms  2Y  -  3, 

(4Y  “  6)/3,  (y  -  7)/2,  and  (Y  +  4)/2*  These  particular  functions 
of  y  are  chosen  because  they  are  easy  to  generate  and  it  is 
possible  to  aoproximate  the  normal  density  very  closely  with  a 
mixture  of  their  densities* 

If  f(y)  is  the  density  of  Y  =  +  U^, 

•  5y^  0  ^7  ^  1 

*5y^  -  1.5(y  -  1)^  1  ^  7  ^  2 

*5y^  -  l*5(y  -  l)^  +  1.5(y  -  zf  2  ^  y  ^  3 
0  y<0ory>3, 
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then  the  densities  of  2Y  -  3,  (4Y  -  6)/3,  (Y  -  7)/2, 

v‘ 

(Y  +  4)/2  are  .5f(.5X  +  1.5),  .75f(.75x  +  1.5),  2f(2x  +  7),  and 
2f(2x  -  4)*  We  use  the  following  mixture  of  these  densities  to 
approximate  the  normal  density* 
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(1)  (2n)“*V*^^  =  .8635[.5f(*5x  +  1.5)]  +  .11506[.75f(.75x  +  1.5)]  + 

.00372[2f(2x  -  4)]  +  .00372[2f(2x  +  7)]. 

Furthermore,  the  expression  on  the  right  is  always  less  them 

_  c  c  2 

(2n)  *^e  *  ;  the  difference  is  plotted  in  Figure  1.  This  means 

that  we  can  generate  X  from  one  of  the  four  densities  on  the 
right  of  (l)  some  98.6  per  cent  of  the  time,  and  the  remaining  1.4 
per  cent  of  the  time  we  generate  X  from  either  the  tail  or  the 
residual  portion  of  the  density.  The  combined  mixture  provides  an 
X  with  density  exactly  (2n)”*^e~*^*  .  The  method  for  handling  X 
from  the  tail  of  the  normal  distribution  is  a  slight  modification 
of  that  described  in  [3],  and  the  rejection  technique  is  used  to 
provide  X  from  the  residual  density  on  -3.5  ^x^  3.5.  See 
Figxire  1.  0\ir  method  for  generating  X  is  then  to  choose  one  of 
six  alternate  ways — from  one  of  the  four  densities  on  the  right  of 
(l),  from  the  residual  density  of  Figure  1,  or  from  the  tail.  The 
details  are  in  this  outline: 

To  generate  a  standard  normal  random  variable  X,  density 
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cp(x)  =  (2n)"*^e"*^*  ,  let  Y  =  ^ 

1.  Wi.th  probability  .8635,  out  X  =  Y  +  Y  -  3 

2.  With  probability  .11506,  £ut  X  =  4Y/3  -  2 

3.  With  probability  .00372,  £ut  X  =  .5Y  -  3.5 
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Figure  1.  G^aph  over  -3.5  <  x  <  3.5  of  the  residual  function 
(2Tt)"*^e'*^^  -.43175f(.5x+1.5)-.086295f(.75x+1.5)-.0074^[f(2x-4)+f(2x+7)]. 

The  area  under  the  curve  ^s  •01353<i7418.  To  generate  a  variate  X  with 
density  the  appropriate  nultiple  of  this  function,  generate  (x,y)  uiiiformly 
from  the  hat-shaped  region  until  (x,y)  lies  under  the  curv^,',  then  put  X~x. 
Choose  (x,y)  uniformly  from  the  hat-shaped  region  by  putting  x  =  7U^  -  3.5, 
y  =  .OOII5U2,  with  probability  .3095558546  or  x  =  ?.9(U^  +  -  1) , 

y  =  .00115  +  .00945U^(1  -  lu^  +  -  l|)  with  probability  .6904441454. 

Efficiency  =  >  or  52  per  cent. 
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4*  With  probability  *00372,  put  X  =  *51  +  2 

5*  With  probability  ,014  “  2  f  9(t)dt  =  •0135347418***  choose  (x,y) 

3.5 

uniformly  from  the  hat-shaped  region  of  Figure  1  until 

Y  <  (2n)"‘^e“*^^  -  .43175f(.5x  +  1.5)  -  .086295f(.75x  +  1,5) 

-  ,00744[f(2x  -  4)  +  f(2x  +  7)] 

then  put  X  =  X, 

CO 

6.  With  probability  2  f  <p(t)dt  =  ,0004652582* ••  generate  pairs 

3.5  . 

X  =  2U^  -  1,  y  =  until  y  <  3,5(12.25  -  2  ln|x|)"*%  then 

put 

1(12.25  -  2  ln|x|)*^  If  x  >  0 

-(12,25  -  2  Lilx|)*^  if  X  <  0  . 

The  values  in  this  outline  are  exact,  so  that  combining  the 
6  steps  with  the  indicated  frequencies  will  produce  an  X  with 
density  exactly  «p(x).  Practically,  of  course,  one  Is  restricted 
by  the  word  length  of  the  computer  for  the  U’s  and  the  probabilities 
In  steps  5  and  6, 


The  fastest  method,  [l],  will  produce  normal  variables  in  the 
IBM  7094  computer  at  the  rate  of  12,500  per  second,  while  the  above 
method  will  produce  8,500  per  second.  In  the  IBM  1620  decimal 
computer  the  fastest  method  yields  50  variates  per  second  versus 
32  for  the  above  method. 
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